Dynamics of waves in ID electron systems: 
Density oscillations driven by population inversion 
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We explore dynamics of a density pulse induced by a local quench in a one-dimensional electron 
system. The spectral curvature leads to an "overturn" (population inversion) of the wave. We show 
that beyond this time the density profile develops strong oscillations with a period much larger 
than the Fermi wave length. The effect is studied first for the case of free fermions by means of 
direct quantum simulations and via semiclassical analysis of the evolution of Wigner function. We 
demonstrate then that the period of oscillations is correctly reproduced by a hydrodynamic theory 
with an appropriate dispersive term. Finally, we explore the effect of different types of electron- 
electron interaction on the phenomenon. We show that sufficiently strong interaction [U(r) 3> 1/ mr 2 
where m is the fermionic mass and r the relevant spatial scale] determines the dominant dispersive 
term in the hydrodynamic equations. Hydrodynamic theory reveals crucial dependence of the density 
evolution on the relative sign of the interaction and the density perturbation. 
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I. INTRODUCTION 

Transport properties of interacting one-dimensional 
(ID) systems keep attracting a great deal of research in- 
terest. Experimental realizations of ID fermionic systems 
include, in particular, carbon nanotubes, semiconductor 
and metallic nanowires, as well as quantum Hall (and 
other topological insulator) edges. Further, ID bosonic 
and fermionic systems can be engineered by using cold 
atomic gases in optical traps of the corresponding geom- 
etry. A standard and powerful theoretical approach to 
interacting ID systems is the bosonization^—. When (i) 
the spectrum is linearized, (ii) backscattering processes 
are neglected, and (iii) the physics near equilibrium is 
explored, the bosonization reduces the original interact- 
ing problem to a Gaussian field theory, thus reducing 
evaluation of physical observables to a straightforward 
calculation of Gaussian integrals. When one (or several) 
of the above three conditions is not fulfilled, the the- 
oretical analysis becomes much more involved. In the 
present paper we will focus on non-equilibrium physics 
of ID fermionic systems in the regime where the spectral 
curvature is of crucial importance. 

Properties of ID interacting systems with spectral non- 
linearity have been addressed in a series of recent theo- 
retical works^— . Here we will consider the time evolu- 
tion of a density pulse created by a local quench in a ID 
fermionic system (that will be assumed to be spinless or 
spin-polarized for simplicity). We will assume that this 
pulse is quasiclassical (i.e., has a characteristic spatial 
extension much larger than the Fermi wave length) and 
sufficiently strong (i.e., contains a large number of elec- 
trons). For not too long times the evolution seems to be 
fully harmless: the pulse splits into left- and right-moving 



parts that separate and move away from each other, ap- 
proximately preserving theirs shape. They key point is 
that the shape would remain strictly unchanged only for 
linear dispersion of excitations, while the non-linearity 
of dispersion leads to a deformation of the pulse. As a 
result, at a certain finite time the pulse tends to "over- 
turn" . The problem to be addressed is what happens 
with the density profile beyond this time. 

The above problem was formulated in Ref. [ll| in the 
context of Calogero model that was argued to describe 
the fractional quantum Hall (FQH) edges (see also a re- 
cent paper Ref. |12) . Using quan tum hydrodynamics ap- 
proach, the authors of Refs. Illlfl2l came to a conclusion 
that a density pulse formed in a FQH edge will evolve 
in a sequence of well separated solitons with a quantized 
charge (equal to v for Laughlin states). 

In the present work we perform a systematic analysis 
of the pulse dynamics for free fermions as well for those 
with different types of interaction. We begin by consider- 
ing a non-interacting case (Sec.[H]). Quantum simulations 
show development of density oscillations at sufficiently 
large times. By analyzing evolution of the Wigner func- 
tion, we show that once the semiclassical phase-space dis- 
tribution overturns (i.e. develops a population inversion 
characterized by three "Fermi momenta"), strong oscilla- 
tions of density are generated in the corresponding region 
of space. The characteristic scale of these oscillations is 
much larger than the Fermi wave length \p. The oscil- 
lations can be understood as Friedel oscillations between 
different Fermi-momentum branches. 

In Sec. IIIII we switch to the bosonization language 
and discuss a connection between free-fermion oscilla- 
tions studied in Sec. |TT] and classical hydrodynamics. 
In the latter class of problems^ oscillating structures 
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are known to develop when shock waves are regularized 
by dispersive terms. We show that although dispersive 
terms arise already within Haldane bosonization formal- 
ism of free fermions with curvature, the dominant terms 
should come from summing the loop expansion. While 
we do not know how to take into account these effects 
systematically, we approximate them by including in the 
classical hydrodynamic equation a term corresponding to 
an upper (for a positive pulse) branch of the particle- hole 
continuum. Solving the corresponding equation (which 
is of Benjamin-Ono type), we show that it yields oscilla- 
tions with correct period (including its spatial variation) 
but with an amplitude several times larger than the right 
one. This shows that the above classical hydrodynamic 
equation does catch some important physics of the devel- 
oping "dispersive shock" of a free-fermion pulse but does 
not represent a fully controllable approximation. 

Section IIVI is devoted to an analysis of the interaction 
effects on the pulse evolution. We consider first the case 
of a short-range interaction and argue that the discovered 
oscillations remain largely preserved, up to two modifica- 
tions: (i) conventional Luttinger-liquid renormalization 
of the Fermi velocity, and (ii) washing out of oscillations 
at long times due to inelastic processes. We turn then 
to the case of a long-range interaction. We show that 
when the interaction decays sufficiently slowly (specifi- 
cally, U(r) 3> 1/mr 2 at large distances r, and m is a elec- 
tronic mass), the leading contribution to the dispersion 
results from the interaction term, and the problem can be 
treated quasiclassically (i.e. loops can be neglected), giv- 
ing rise to a classical hydrodynamic equation. The evolu- 
tion of the pulse according to such an equation depends 
crucially on the sign of the pulse and the sign of the in- 
teraction (or, more precisely, on the relative sign between 
them) . When the interaction is repulsive and the density 
pulse is downward, oscillations develop similarly to the 
case of free fermions (or short-range interaction). The 
period of oscillations gets however parametrically larger. 
On the other hand, for an upward pulse (and still assum- 
ing a repulsive long-range interaction), the pulse splits in 
a sequence of "solitons" (whose charge is in general not 
quantized, except for the case of 1/r 2 interaction). 

Section [V] contains a summary of our results and a 
discussion of prospective research directions. 



II. FREE FERMIONS 

In this section and in Sec. IIIII we study the evolution 
of a "quasiclassical" density disturbance of the Fermi sea 
of free fermions. We begin (Sec. Ill A|) by formulating the 
problem and performing its numerical modelling which 
shows emergence of density oscillations after the time 
corresponding to overturning of the initial packet. In 
Sec. Ill Bl we solve this problem analytically by using the 
straightforward ( "fermionic" ) approach. Specifically, we 
demonstrate that, once the dispersion induces a popula- 
tion inversion within the pulse, phase-space oscillation of 
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FIG. 1: Setup. Density disturbance in the Fermi sea is created 
by the application of a potential U(x) which is then switched 
off at t = 0. 

the Wigner function give rise to density oscillations. Ana- 
lyzing the resulting density oscillations, we find a perfect 
agreement with the results of numerical simulations of 
Sec. Ill Al Finally, in Sec. IIIII we make a link to "dispersive 
shocks" in the classical hydrodynamics. We show that 
when the dispersive term corresponding to the appropri- 
ate branch of the particle-hole spectrum is incorporated 
into classical hydrodynamic equations, the latter repro- 
duce correctly the period of emerging oscillations (but 
considerably overestimates their amplitude). 



A. Formulation of the problem and numerical 
simulations 

The problem that we address is formulated in a rather 
simple way. We assume that a non-uniform fermionic 
density was created by application of a smooth (on the 
scale of Af) and relatively weak (compared to the Fermi 
energy ep) external potential U{x) to the unperturbed 
Fermi sea (Fig. [TJ. The system at t < is in its ground 
state characterized by the fermionic density 



Po ( x ) = t£U. = - y / p 2 2 mU{x). (1) 

Here, is the Fermi momentum at infinity. Note that 
all the corrections to the semi-classical result ([T]) are ex- 
ponentially small as long as U (x) is smooth on the scale 
of Af- For transparency of discussion we assume that 
the density pulse has a shape of a single hump (as shown 
schematically in Fig. [T]), i.e. that U(x) < and has a 
single minimum at x = 0. 

At t = the potential is suddenly switched off, which 
results in the appearance of a non-equilibrium state and 
subsequent propagation of the density perturbation cre- 
ated by U(x). Our goal will be to explore this den- 
sity evolution at sufficiently long times. We will assume 
that the number of particles within the initial pulse is 
large, Ax Ap ^> I, where Ax and Ap are the char- 
acteristic extension of the pulse and its amplitude, re- 
spectively. The interesting physics will emerge at times 
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FIG. 2: Evolution of the fermionic density (measured in fermions per site) as obtained from direct numerical simulations of 
free fermions on a lattice. The unit of length is set by the lattice constant. Left: The initial density pulse has Gaussian shape 
with the amplitude Ap ~ 0.05 and the dispersion Aa; ~ 102.6, comprising 12.8 electrons. Middle: Density profile at t slightly 
smaller than the "shock" time t c . Only right-moving part is shown. The front edge of the pulse has become steep but the 
overturn has not yet occurred. Right: Fermionic density after the "shock", t « 5t c . Density ripples develop at front edge. 



t > t c ~ mAx/Ap, when the semiclassical phase-space 
distribution overturns. 

Since the fermions are free and their state was prepared 
in the coherent manner described above, the full informa- 
tion on the quantum state of the system is encoded in the 
Wigner function 



f(X,p;t) = / dye-^(i/>+(X-it)il>(X + ;t)) 



(2) 
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satisfying at t > the (exact!) Boltzmann equation 

d t f(X,p;t)+pd x f(X,p;t) = 0. (3) 

In this equation we have set the particle mass to unity. 
Let us note that the mass m will enter our results only 
through the overall time scale. Thus, dependence on m 
can be eliminated completely by measuring t in units 
of t c . 

Equation (E| is trivially solved, yielding 

f(X,p;t) = f(X-pt,p;t = 0) = f (X-pt,p) . (4) 

Therefore, once we know the Wigner function of the ini- 
tial (t = 0) state, the Wigner function of the evolved 
(t > 0) state can be immediately obtained. 

Let us examine now the Wigner function of the ini- 
tial state. Semiclassically one would expect that fo(X,p) 
takes value unity for the occupied electronic states 
that are below the position-dependent Fermi momentum 
Pf{X) and is zero for empty ones (above pp(X)): 



fo(x,p) = e(p F (x)- P 2 ) 



(5) 



In this approximation, the time-dependent state of the 
fermions after the quench (i.e. at t > 0) is fully char- 
acterized by the Fermi surface pp(x,t) separating occu- 
pied and unoccupied single-particle states in the phase 
space and satisfying the Euler equation (we concentrate 
on the Fermi surface for the right-moving particles with 
p F (x) > 0) 



While capturing correctly the physics at small times, 
the Euler equation © suffers from the shock-wave phe- 
nomenon. Specifically, for arbitrarily smooth initial con- 
ditions, the curvature of the electronic dispersion relation 
e (p) — P 2 /2, makes the Fermi surface pf{x) multival- 
ued at large enough times (t > t c ) and leads to the ap- 
pearance of infinite spatial gradients of fermionic density 
(Fig. El middle and bottom panels). This suggests that 
the simple semiclassical description ©, © may become 
insufficient beyond the time t c when the shock occurs, 
raising the question of what happens with the density 
profile at t > t c . 

We have performed direct quantum simulations of this 
problem by using a tight-binding free-fermion model. 
Figure [5] demonstrates the density evolution from initial 
state at t = (left panel) to state at certain time after 
the shock (approximately five times larger than t c ). At 
t > only the right-moving part of the density pulse 
is shown. A full movie of density evolution is available 
online 14 . The initial density perturbation was Gaussian 
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(7) 



d t pF(x,t) +p F (x,t)d x p F (x,t) = 0. 



(6) 



with dispersion a of about 100 lattice sites and contained 
N s» 12.8 electrons. The density of the underlying Fermi 
sea is 0.2 fermion per site, so that the cosine-shaped dis- 
persion relation of the tight-binding model can be well 
approximated by a parabola. 

For convenience of the reader the snapshot of the den- 
sity at f ~ ht c is also shown in the top panel of Fig. El 
Comparing the exact quantum result (Fig(3l top) to the 
naive semiclassical result dictated by the Euler equation 
(Fig. El bottom), we see that the shock gets regularized 
via the onset of pronounced density oscillations at the 
front edge of the pulse. It is important to emphasize that 
the period of those oscillations is controlled by the am- 
plitude Ap of the density perturbation and is thus much 
larger than Xp. (We will perform a detailed quantita- 
tive analysis of the oscillation period below in Sec. Ill Bl ) 
From this point of view, the developing oscillations may 
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FIG. 3: Top: Snapshot of fermionic density t w 5t c as ob- 
tained from direct numerical simulations of free fermions on a 
lattice. The initial density pulse was Gaussian (see main text 
and Fig. [2]). The dots show the semiclassical result derived 
in Sec. Ill Bl Middle: Quasiclassical phase-space distribution 
as obtained from the Euler equation at t m 5i c . In between 
trailing and leading edges xt and xi the Fermi surface is mul- 
tivalued with branches , i = 1, 2, 3. Bottom: Naive quasi- 
classical approximation to the density profile as obtained by 
momentum integration of the phase-space distribution shown 
in the top panel. 



be considered as quasiclassical: their characteristic scale 
is much larger than \p- Thus, a smooth initial density 
stays smooth at scale Xf also at times after the "shock" . 

We thus face an apparent contradiction: the density 
profile remains "quasiclassical" (smooth on the scale of 
Xf) after the shock but develops strong oscillations that 



are not caught by the quasiclassical approximation based 
on Eqs. ([5]), ©. The resolution of this "paradox" is re- 
lated to the fact that Eq. §5§ is not the fully correct 
semiclassical (in the above sense) approximation for the 
Wigner function of fermions in a smooth potential well. 
As was pointed out in Ref. flH instead of having abrupt 
drop from 1 to at Fermi momentum, fo(X,p) as a func- 
tion of p develops oscillations near pp{X). Those oscil- 
lations can be considered as a semiclassical effect in the 
sense that their form knows nothing about A p and is con- 
trolled solely by the derivatives of pp(X). In Sec. Ill Bl we 
give a detailed account of the oscillations in the Wigner 
function and of their implications for the density evolu- 
tion. 



B. Wigner function of fermions in a potential well 
and density oscillations 

The Wigner function fo(p,X) of the initial state sat- 
isfies the equation 

d x d y f (X, y) - [U{X + ) - U(X_)] f (X, y)=0, (8) 

where the coordinate y is conjugate (in the sense of 
Fourier transformation) to the momentum p, cf . Eq. ^ . 
As we are interested in the behavior of fo(p,X) close to 
one of the Fermi edges p = =Lpoo we can replace d y by ipoo 
(we concentrate here on the right Fermi edge). This cor- 
responds to taking the limit p^ — » oo while keeping the 
profile Pf(x) — p^ fixed. Solving the resulting equation 

dxfo(X, y) - i {pf(X+) - p F {X-)] f (X, y) = , (9) 

with the condition at infinity 

fo(X = -oo, y) = J |V^9 ( Poo - p) (10) 

and transforming the result to momentum space, we find 
in agreement with Ref. [TH 



fo(X,p) 



dy 



-iS[y;X,p] 



2iti(y — iO) 
S[y;X,p]=py- dX'p F (X') 



(11) 

(12) 



Note that our approach is slightly different from that of 
Ref. [ToT : we consider an equilibrium state in a potential 
U{x), while the authors of Ref. UH construct a coherent 
state by acting on the homogeneous Fermi vacuum with 
an exponential of a bilinear in fermionic operators. 

The general structure of the Wigner function can be 
inferred from Eqs. pip . (|12l) by performing the integra- 
tion with making use of the saddle-point method. The 
saddle-point equation for the action (fT"2")) 
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FIG. 4: Initial-state Wigner function fo(X = 0,p) as ob- 
tained via numerical integration of The same Gaussian 
density profile as used in the quantum simulations presented 
in Sec. Ill Al and in Fig.[2]was assumed. The observed six oscil- 
lations correspond to 6.4 right-moving particles in the pulse. 
Vertical lines in Fig. [4] mark the momenta satisfying the con- 
dition for oscillation maxima, S[X,p] = — 3n/4 mod 2n. 
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has no real- valued solutions for p < p^x, or 

PF (X + y/2)+p F (X-y/2) 



P > PmaxO) 



max • 

y 



(14) 



In these parts of the phase space the integral (fTTj) is con- 
trolled by the singularity at y = 0, and fa(X,p) can be 
well approximated by Heaviside 0-function. 

On the contrary, for p^ < p < p max (X) there ex- 
ist (at least) two solutions y = ±y*(X,p) to the saddle 
point equation providing an oscillatory contribution to 
the Wigner function 



5/o(X,p)cx5R 



Ae 



-iS[X,p] 



(15) 



Here A is the coefficient controlled by the fluctuations 
around the saddle points. The phase of oscillations is 
given by S[X,p] — S[y*(X,p);X,p\. It defines the pe- 
riod of oscillations with the momentum via S[X, p+Ap] — 
S[X,p] — 2n. Counting the powers of imaginary unit in 
the saddle point integration one easily finds that the max- 
ima of the Wigner function appear at S[X,p] = —3ir/4 
mod 2tt. 

In a sufficiently close vicinity of the local Fermi surface, 
\p — pf{x)\ <C Pf(0) — Poo one can locally approximate 
Pf(x) by a parabola and express the Wigner function 
fo(x 7 p) in terms of the Airy function^. This approx- 
imation is however insufficient for our purposes, as the 
density oscillations will originate from Wigner function 
oscillations at all scales \p — pf(x)\ ~ Pf(0) — Poo. 

The oscillatory behavior of f (X,p) is illustrated in 
Fig. [4] where we plot f(X = 0,p) (calculated via numer- 
ical integration of pip ) as a function of momentum. To 
generate the plot we have assumed the Gaussian density 
used in the quantum simulations presented in Sec. Ill Al 



FIG. 5: Phase-space oscillations of the Wigner function 
fo(x,p) for the initial state (the same Gaussian density hump 
as in the previous figures; the corresponding p F (x) is shown 
by dotted line). Blue solid lines are contours of constant ac- 
tion (which determines the phase of the Wigner function), 
S[X,p] — — 37r/4 mod 2-7T. Dashed line represents the border 
p max (x) of the region of developed oscillations. 



and in Fig. [2] A straightforward analysis of the saddle- 
point equation shows that, upon variation of momentum 
from p = p m ax(0) = p F {0) to poo, S[X = 0,p] varies 
monotonically from to —2ttN, where 



N 



1 

2^ 



dx(p F {x) - Poo) 



(16) 



is (generally non-integer) number of particles in the pulse. 
Accordingly, the Wigner function of Fig. |4] shows 6 os- 
cillations corresponding to approximately 6 right-moving 
particles. Vertical lines in Fig. U mark the momenta sat- 
isfying the condition S^Ajp] = — 37r/4 mod 2ir. 

The overall behavior of fo(X,p) in the phase space can 
be conveniently represented by lines of constant action 
^[V.p]. For the case of Gaussian density the correspond- 
ing pictures is shown on Fig. \5\ The green dashed line 
here represents the border p max {x) of the region of devel- 
oped oscillations. The solid blue lines are the lines of con- 
stant action S^V,;?] = —in/A mod 2ir. Finally, the dot- 
ted line shows the x-dependent Fermi level. The "topol- 
ogy" of the plot can be understood on general grounds 
and does not depend on the specific density po(x). At 
large X the action is a steep function of momentum p 
and, when taken modulo 2-7T, acquires any given value 
many times. Among the lines of constant action coming 
from X — +oo, exactly [N] (this denotes the integer part 
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of N) lines cross the p-axes and flow to X = — oo, while 
other lines end up on the border p = p max (x) where the 
solution y*{X,p) to the saddle point equation becomes 
complex. 

Let us now discuss the implications of the above results 
for the fermionic density which is equal to the integral 
over momentum of the Wigner function. In the initial 
state the density is insensitive to the oscillations of the 
Wigner function. Indeed, one can observe that the con- 
tour of momentum integration (vertical line on Fig. [5]) 
crosses many contours of constant S[X,p] but does not 
touch any of them, so that there is no stationary-point 
contribution to the integral. In fact, evaluating the p in- 
tegral of Eq. (fTTj). we find that in the considered (large 
Pf) approximation the equality po(x) = pf{x)/it is ex- 
act. 

The evolution of each contour line of the action is gov- 
erned by the Euler equation ©. The behavior of the 
Wigner function after the shock, t > t c , is illustrated by 
Fig. [6l Now the vertical lines do touch the contours of 
constant action. A touching point becomes the saddle 
point for the integration over momentum and the oscilla- 
tions in f (X,p) start to contribute to the density. This 
implies a maximum in the density each time the inte- 
gration line touches the contour line corresponding to 
S[X,p] = -3tt/4 mod 2tt. 

The brown curve in the upper part of Fig. [5] shows the 
density obtained via numerical integration of the Wigner 
function This is almost indistinguishable from the 

result of first-principle quantum simulations (top panel of 
Fig. 13]). We observe that the positions of the maxima of 
the density are in accord with the above argument based 
on the saddle-point approximation. 

One can now determine the period of the density mod- 
ulations. Let us focus on the region closer to the front 
edge of the pulse where the oscillations are clearly gov- 
erned by a single harmonics, see Fig. [3J and are per- 
fectly described by the saddle-point argument as shown 
in Fig. [51 As this figure further illustrates, the integration 
contour touches the contour lines near p = p max (A) and 
the period of the density oscillations is set by S[X,p = 
Pmax(A)]. In this regime we have 

y*(X,p max (X))~2X, 

p^iX)***™^. (17) 

Thus, S[X,p ma _x(X)] ~ (pf(0) — Poc)X and we immedi- 
ately infer the period 



Pf(0) - Poo 

It is not difficult to generalize this argument to the 
region closer to the top of the pulse. Assuming for sim- 
plicity that the time that has passed after shock is of the 
order of the shock time t c ~ mAx/Ap, we find that the 
characteristic spatial scale for the first few oscillations 
(just to the right of the maximum of the pulse) is larger 




3000 3500 4000 

FIG. 6: Top: Density profile after the shock as obtained 
by numerical numerical integration of the Wigner function 
(|lip . The result is almost indistinguishable from that of first- 
principle quantum simulations shown in Fig. [3] Bottom: The 
same as in Fig. [5] but after evolution beyond the shock. All 
parameters, including the evolution time, are the same as in 
the upper plot and in Fig. [3] The vertical lines touching the 
constant-action contours mark stationary-point contributions 
to the momentum integral of the Wigner function and thus 
determine maxima of the density oscillations in the upper 
plot. 



than (|18p by a factor ~ N 1 / 3 (for the parameters used 
in our plots this is approximately two). Some compli- 
cation comes from the fact that in this region a super- 
position of oscillations originating from different regions 
in phase space takes place. Indeed, as is clearly seen in 
the lower panel of Fig. [BJ constant-action contours that 
do not terminate at p niax (x) may have two points with 
infinite slope, and each of them will give a stationary- 
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phase contribution when the momentum integration is 
performed. This superposition explains a somewhat ir- 
regular oscillation pattern in the corresponding spatial 
region. 

The obtained oscillations can be interpreted as Friedel- 
type oscillations between different branches of the 
Fermi momentum (that becomes multivalued after the 
"shock" ). In particular, in the front region, the upper two 
branches are close to the maximum value, po(0), while 
the lower branch is essentially equal to p^ , which yields 
exactly Eq. (flgj). 

III. HYDRODYNAMICS OF FREE FERMIONS 

The analysis of the previous section provides a detailed 
description of the evolution of coherent perturbation in 
the density of free fermions. The analysis is complete 
and, as was also confirmed by numerical simulations, es- 
sentially exact under our basic assumptions. It is appeal- 
ing however to try to formulate a hydrodynamic descrip- 
tion of the problem which, in contrast to the fermionic 
approach utilizing the notion of Wigner function, would 
involve as fundamental objects only the density and the 
velocity of the electronic fluid. Indeed, the hydrodynam- 
ics (bosonization) constitutes a convenient and powerful 
framework for the discussion of interaction effects (to be 
considered in Sec. IIVI below) which are otherwise hard to 
access. 

Usually, hydrodynamics rests on the assumption of lo- 
cal equilibrium forced by the particle collisions which 
wash out any features in the particle distribution func- 
tion. In the present problem no such equilibrium ex- 
ists. Moreover, we saw above that the oscillating behav- 
ior of f(p) is crucial for the density ripples observed in 
the shock region. Hence, one can expect that the evo- 
lution of the quantum-coherent many-particle state can 
not be controllably described by classical equations of 
hydrodynamic type. Despite this fact, one can ask if 
it is possible to design phenomenological hydrodynamic 
equations which would capture qualitative features of 
the true density evolution. In this section we show that 
this is indeed possible and such phenomenological equa- 
tions provide important insight into the physics of non- 
equilibrium many-particle system. 

In our search for hydrodynamics it is convenient to 
start from the Euler equation ([6]) corresponding to the 
neglect of all oscillatory features in the Wigner func- 
tion f Q (X,p). Combined with analogous equation for 
the Fermi surface of left electrons at t < t c it can be 
rephrased in terms of the mean density and velocity of 
the fluid as 

(2 2 2 \ 

y + 1 2~) =0 ' (19) 

p{x) = J T^ffrP) > «(*) = ^) / ^ Pf( ~ x ^ -( 2 °) 
Of course, equations (fH?|) suffer from the same shock 



phenomenon as the original equation ([6]). Phenomeno- 
logically, we would like to add some terms to Eq. (fT9")l 
regularizing the shock instability. We know that reg- 
ularization goes trough the onset of density ripples in 
the shock region. This phenomenon is well known in 
hydrodynamics and usually referred to as "dispersive 
regularization"—. It takes place when the shock caused 
by non-linearity gets regularized by higher order deriva- 
tives consistent with the time reversal invariance of the 
equations. The classical examples are the Korteweg-de- 
Vries (KdV) and Gross-Pitaevskii equations. It is the 
requirement of time reversal that makes the "dispersive 
regularization" very different from "dissipative regular- 
ization" achieved by the introduction into the system of 
some type of viscosity. Examples of the latter type are 
Navier-Stokes and Burgers equations. 

Let us now point out that Eqs. (1191) appear in the the- 
ory of Fermi gas in yet another context and with slightly 
different meaning. Specifically, the standard bosoniza- 
tion procedure applied to the fermions with quadratic 
spectrum leads to Hamiltonian— ~— 

Here, p and v are operators with the commutation rela- 
tions 

[p(x),v(y)]=-i6'(x-y). (22) 

The Hamiltonian and the commutation relations imply 
the operator equations of motion usually referred to as 
"quantum Euler equations" 

/ ^2 2 -~2 \ 

d t p+d x (pv)=0, d t v + dj^- + =0. (23) 

One can now see that there exist two sources of correc- 
tions to hydrodynamic equations (|19p. First, there can 
be corrections to the quantum Hamiltonian (|21[) missed 
by the bosonization in its simplified form. If present, they 
would yield a direct contribution to the quantum Euler 
equations (|23[) . Second, passing from quantum equations 
([23]) to identically-looking classical equations (p~9|) implies 
averaging of the former over the quantum state. In the 
functional integral formulation of the problem, classical 
equations (|19l) correspond to the saddle-point treatment 
of the functional integration. However, loop corrections 
can also contribute to the average density and current 
and generate new terms in (|19p . Below we discuss both 
aforementioned effects. 



A. Correction to Hamiltonian 

Let us first explore corrections to the Hamiltonian (|21[) . 
For a while we put the loop corrections aside (we will 
return to them in Sec. MI B| and thus make no distinction 
between classical and quantum equations of motion. 



8 



We start with the Haldane's theory^ that accounts 
for a discrete nature of particles as well as for their 
spectrum? 1 . Within this model, the fermionic operator 
is represented by an infinite sum 



*(x) 



' pe 



i6(x) 



— iln4>(x) 



(24) 



where the bosonic fields have the standard commutation 
relations 



[0(30,000] = ~sgn(x-x'), (25) 
and are related to the velocity and density fields as 

v(x) = dj{x) , p{x) = dj(x) . (26) 
After substituting Eq. (f24|) into the free Hamiltonian 



H 



dx^ + (x)V 2 ^(x) , 



one obtains 
H = 



J l odd 



1P| 

4 p 



(27) 



(28) 



This result contains an infinite summation over odd inte- 
gers I, and formally diverges. To properly define this se- 
ries, one needs to regularize the divergent sums. This can 
be achieved by describing the series as an expansion of an 
analytic function of some argument (y(z) = yiz ). A 
series of this type has to be summed within the range 
of its convergence and then analytically continued to 
z = 1 . Bearing such a procedure in mind and comparing 
Eq. (j2"5)l with Eq. (|2"T1) , we establish that such regulariza- 
tion implies 



^2 1 = 1 and 12 = V 3 - 



l odd 



Thus, one obtains 



H = j dx[ -ptf 



l odd 



6 8 p 



-2 



(29) 



(30) 



The first two terms in Eq. (|30p arc contained in the 
Hamiltonian (|28[), However, the last term in Eq. (|30|) 
represents the gradient corrections that are beyond the 
Hamiltonian (|2"8")l . Note that though the derivation out- 
lined above may not appear rigorous, there is no ambi- 
guity in determining the last term in Eq. (|30l) . Indeed, 
the regularization procedure we use is fully determined 
by the first two terms in the Hamiltonian. Therefore, the 
coefficient in front of the third term is unambiguously de- 
termined. The Hamiltonian ([30]) leads to the following 
equations of motion 

d t p + d x (pv) = 0, d t v + vd x v + d x w = . (31) 



P(x) 
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FIG. 7: Results of numerical solution of hydrodynamic equa- 
tions (|31[l . (I32|l . The initial density pulse was gaussian with 
the parameters used previously in Sec. ITT] The period of os- 
cillations induced by the shock is determined by t he comp eti- 
tion between dispersion and non-linearity, Sx ~ ^/Xp/Ap. It 
is parametrically smaller then the one found in Sec. HTl within 
a direct analysis of the free-fermion problem. 



Here 



(32) 



7T 2 /5 2 1 2 „ 1 2 

w = — -d x \ogp- - (d x log p) 



is the enthalpy of the Fermi gas. The first term in Eq. ([32|) 
is the pressure of a homogeneous Fermi gas, while the last 
two terms describe the cyclotronic pressure that accounts 
for the finite density gradient. Interestingly enough, this 
latter contribution is quite universal and appears also in 
the Madelung fluid 1 as well as in the hydrodynamic form 
of the Gross-Pitaevskii equation^. The presence of the 
finite gradient terms stabilizes the classical equation of 
motions. Thus, it is natural to ask a question whether 
Eqs. (|3"TT) . (|32l) arc sufficient to describe the evolution of 
the density pulse discussed in Sec. Ill Al 

To answer this question, we simulate the evolution of 
the density pulse ([7]) in accordance with Eq. (|32|) . The 
results of this analysis shown in Fig. [7] clearly indicate the 
formation of a region of oscillations in the density profile. 
However, the period of the oscillations is parametrically 
different from that obtained from the direct quantum- 
mechanical solution of the free-fermion problem of Sec, fill 
see Appendix [S] for details. Indeed, the spatial scale of 
the oscillations is determined by the competition between 
the non-linearity and the dispersion. In the present case, 
a simple estimate gives 



Sx 



1 /Af 

y/paoAp V Ap 



(33) 



This is smaller by a factor (XpAp) 1 ' 2 -C I than the re- 
sult (fT8|) of the direct solution of the free-fermion prob- 
lem. We thus conclude that equations (I5T|) . (|32p yield 
a parametrically wrong scale for the density ripples: the 
dispersive term in Eqs. (I5D1 . (|32p is too weak. Thus, 
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in our search for hydrodynamics we must resort to loop 
corrections to the equations of motion. 

Before passing to the analysis of loop corrections, let us 
make the following comment. While the dispersive term 
in Eq. pop turns out to be parametrically small in com- 
parison to quantum effects for free fermions, such a term 
(with a parametrically enhanced prefactor) will become 
a dominant dispersive term for the case of electrons with 
finite range interaction with a sufficiently large interac- 
tion radius, see Sec. IIV Al Consequently, the semiclassi- 
cal analysis of Eq. pop (with an appropriately modified 
coefficient of the last term) performed above and in Ap- 
pendix [X] will become a controllable description in that 
case, as discussed in Sec. IIV Al 
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B. Loop corrections 

Let us follow our phenomenological approach and try 
to guess the form of the loop corrections to the enthalpy 
(|32|) on the basis of our knowledge of the characteristic 
scale of the ripples. It is easy to see that to produce 
the correct period of the density oscillations the correc- 
tion should scale as first power of momentum. A sim- 
ple term of the form d x p is not acceptable as it would 
break the symmetry with respect to the spatial inver- 
sion. This symmetry can be saved however by inclusion 
of the Hilbert transform H, 

Sw ~ Hd x p . (34) 

By definition, in momentum domain the Hilbert trans- 
form acts according to Hpk — —iirsign(k)pk- From now 
on we will reserve a special notation Ai for the operator 
Hd x . The reason for such a notation will become clear 
in Sec. IIVI In momentum space 

Mpk = n\ k \Pk ■ (35) 

The enthalpy corrections of the form (|34|) were first 
suggested by Jevicki^ 2 - in his study of the p 3 -theory de- 
fined by the the Hamiltonian ([2~Tj) . He pointed out that 
such a theory contains two single-particle branches. In 
fermionic language these correspond to the electron and 
hole parts of the spectrum, 

k 2 

e p (h) = k F\k\ ± — , (36) 

with subscripts p and h referring to particles and holes, 
respectively. Further, it was observed in Ref. that 
each of the Lagrangians 

* = 3-y (38) 

when treated semiclassically (i.e. on the saddle-point 
level), reproduces correctly the dispersion relation for 



FIG. 8: Results of numerical solution of the semiclassical hy- 
drodynamic equations (|31[) . (|39|) . Vertical lines mark the po- 
sitions of density maxima in the exact quantum-mechanical 
result for free fermions, see Sec. HTl 



the corresponding branch of excitations (j3"r?)) . In other 
words, Eq. (|3"Tf is an effective semiclassical theory that 
takes into account explicitly quantum corrections of the 
original cubic theory (|21l) . 

The term with the Hilbert transform in Eq. (f3"Tf gives 
rise to k 2 correction to the linear spectrum of the conven- 
tional bosonization, see Eq. (131)1) . It is known that loops 
in the perturbative diagrammatic treatment (in the con- 
text of equilibrium problems) of the Hamiltonian (|2I[) 
lead indeed to such an effec t 23 i 24 . Specifically, with loops 
taken into account, the support of the bosonic spectral 
weight in the (u>, fc)-plain, which is just a line ui = kpk 
for the linear electronic spectrum, starts to receive a fi- 
nite width 25 of order k 2 . We thus see, that the inclusion 
of the term (|34|) into the enthalpy is a natural way to 
simulate the effect of loop corrections. 

Motivated by this findings, we try to apply the effec- 
tive semiclassical Langrangians (|37|) to our problem. At 
this point, a question naturally arises: which of the two 
Lagrangians L p , Lh should we choose [i.e., which sign 
should we choose in Eq. (|37|) ]? We argue here in the 
following way. Let us assume that the original density 
perturbation is positive, i.e, has a form of a hump as 
shown in Fig. [TJ Such a perturbation can be obtained by 
generating particle excitations on top of a homogeneous 
vacuum state. Therefore, we choose the Lagrangian L p 
as appropriate in this situation. Similarly, in the case 
of a dip-like (i.e, negative) density perturbation, the La- 
grangian Lh should be taken. This choice is by no means 
innocent, as will be discussed in more detail below. 

For definiteness, we consider a hump-like excitation 
(as was also done in Sec. |llj and thus the Lagrangian 
L p . Corresponding hydrodynamic equations (j3Tj) with 
the enthalpy 

w = - \dl \og P - 1 - (d x log pf+Ap (39) 
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are of the Benjamin-Ono type. They were studied previ- 
ously in the literature in the context of Calogero model 
(see Sec. HVBj) . In Sec. MI Cl we analyze these equations 
and compare the outcome to the result of the fermionic 
solution presented in Sec. Ill Bl 



C. Non-local hydrodynamics of free fermions 

As a first step of our analysis of the classical hydro- 
dynamics defined by Eqs. (j3"Tj) . (|39|) . we have performed 
their numerical simulations for the initial density used 
previously in our quantum computations. The result is 
shown in Fig. [5] The dashed vertical lines mark the po- 
sitions of the maxima in the exact fermionic density ana- 
lyzed in Sec. UU (the same lines as in Fig. [6]). We observe 
a very good agreement between the free-fermion problem 
and the classical hydrodynamics I3T) , (1551) in the period 
of the oscillations induced by the shock. This agreement 
becomes essentially perfect close to the front edge of the 
impulse. 

To explore analytically the oscillations emerging in 
the hydrodynamics Eqs. (j3"Tj) . (|39)l at times exceeding 
the shock time t c , we employ the Whitham modulation 
theoryi 3 .. Within this approach, one considers the so- 
lution to hydrodynamic equations in the shock region 
(the interval between the points xi and x t of Fig. [3] as 
a periodic single-phase wave with slowly modulated pa- 
rameters (wave vector, frequency amplitude, etc.). The 
modulation equations for those parameters are obtained 
from the Lagrangian averaged over a period of oscilla- 
tions. For the Lagran gian L p the single-phase periodic 
wave was found in Ref. I36l: 



[x,t) = p x - -ft + $(0) , 



kx-Lut. (40) 



Here po and 7 represent mean density and the current 
in the wave; k and ui are the wave vector and frequency, 
and 27r-periodic function $(#) is defined by its derivative 



tanh a 



■ If sinha 
$ = — 1 

2tt \ cosh a — cos ( 
47rfc 3 pf 

k 4 pj + 4ir 2 k 2 pf - 4(/c 7 - upo) 2 ' 



(41) 
(42) 



Pi = Po + 7T ■ ( 43 ) 



The density in the wave p = d x 4> reads 
k sinh a 



P = Pi 



2n cosh a — cos ( 



(44) 



The parameter a controls the amplitude of the periodic 
wave 



A = 



A 



max Pmxn 



27rsinha 1 



(45) 



as well as its shape. In the limit a> 1 Eq. (|4"4")l reduces 
to weak harmonic oscillations, while in the opposite limit, 



a <C 1, one gets a train of well separated solitons, each 
of them carrying exactly one electron. 

The modulation equations for parameters 7, po, k and 
lo were derived in Ref. [26| where a specific (Lorentzian 
form of the original pulse was used (see also Refs. I27ll2 
for the discussion of modulation equations for a closely 
related Benjamin-Ono equation). In Appendix [B] wc 
present a general derivation of modulation equations. 
The result is conveniently presented in terms of four Rie- 
mann invariants it, satisfying 



d t Ui + UidxUt = . 



0. 



..4. 



The parameters of the wave are given by 



«2 - Ul 



PO 

1 = 



CO 

u 3 - 



k 

= \ (ui - uf) 

U2 + Ul - Uq 



2tt 



4?r 



(46) 

(47) 
(48) 

(49) 
(50) 



The modulation equations should be supplemented by 
the boundary conditions at the ends of the shock region 
x t and xi. These conditions consist of the requirement 
that the average density and current po and 7 match 
those dictated by the Euler equation © in the regions 
without population inversion (x < xi and x > Xt). We 
thus have 



(3) , 
P F +Pao 

Pa = oZ — ' 7 = 

(!) 1 
P F +P00 

P0 = ^r- ' 7 



pf 



pL 



(#) 



47T 
2 



x = x t , (51) 



4tt 



where we used the notation Pp introduced for three 
branches of the Fermi momentum in Sec. [Til see Fig. [3] 
The solution of the equations for Riemann invariants with 
these boundary conditions is given by 



, x = xi , (52) 



(0 

Ui = py , 



* = 1,2,3, 

"0 = -Poo 



(53) 
(54) 



The modulation theory described above reveals a deep 
connection between the hydrodynamic system ([3T1 [39]) 
and the free fermions. Indeed, according to Eq. (|53l) . 
the Riemann invariants m, i = 1,2,3, characterizing hy- 
drodynamic density oscillations in the shock region are 

(i) 

exactly equal to three branches p F of the Fermi surface 
of free fermions in the population-inversion regime. Fur- 
thermore, the equipotential lines of the action S[X,p] 
that played a central role in our "fermionic analysis" 
of the density ripples evolve according to exactly the 
same Euler equation © as that for Riemann invariants, 
Eq. 06]). 
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Equations (|47l) . ([53)1 allow us to make a precise state- 
ment on the period of oscillations: 



6 x 



2n 



(2) 

P F 



(55) 



Thus, we see that close to the front end of the pulse the 
period is 



Sx = — — . 
Ap 



(56) 



and coincides exactly with that of the density oscillations 
for the exact quantum-mechanical solution of the free- 
fermion problem. It is easy to see that near the top of the 
pulse the period is larger by a factor ~ N 1 ^ 3 (assuming 
for simplicity that the time t—t c that has passed after the 
shock is of order t c ), again in agreement with the analy- 
sis of Sec. Ill Bl Equation (|55|) is fully consistent with the 
interpretation of the oscillatory structure as Friedel os- 
cillations between different Fermi-momentum branches. 

The following comment is in order here. As has been 
explained above, when the hydrodynamic theory (1311) . 
([55)1 is used to describe the behavior of free fermions, the 
choice of the particle branch L p in Eq. (|37l) captures es- 
sential features of the evolution of a density hump, while 
the hole-branch Lagrangian Lh is appropriate for a den- 
sity dip. On the other hand, if we would try to apply, e.g., 
Lh for a density hump, it would fail completely. Specifi- 
cally, it would predict the formation of the solitonic train 
in front of the running pulse, i.e. a decomposition of the 
initial density perturbation into well separated solitons 
(cf. Sec. HVBp . which never happens for free fermions. 
It remains an open question whether there exists an im- 
proved hydrodynamic theory that describes evolution of 
the free-fermion density perturbation consisting of a com- 
bination of humps and dips. In this context, it is also 
worth reminding the reader about the following. While 
the classical hydrodynamics analyzed in this section per- 
fectly reproduces the period of free-fermion oscillations 
induced by a shock, it considerably overestimates their 
amplitude. If there exists a better hydrodynamic descrip- 
tion of this problem, one might hope that it would be free 
also of this drawback. 



IV. INTERACTION EFFECTS 

In the previous sections we discussed the evolution of 
the density perturbation in the free electron gas within (i) 
the exact "fermionic" approach and (ii) the phenomeno- 
logical hydrodynamics. We concluded that, in the latter 
formalism, quantum loop corrections are crucial in de- 
termining the character of the dispersive regularization 
of the shock. They can be modeled qualitatively by the 
non-local term (|34l) in the enthalpy of free fermions. At 
this level, the bosonized Hamiltonian for free fermions 
(that corresponds to the effective Lagrangian 1571) takes 



the form (different for particle- and hole-like perturba- 
tions) 



H 



P (h) 



dx 



-pv 



ir 2 p 3 
6 



8 p 



1 - 



(57) 



The term with p\ coming from Haldane bosonization pre- 
scription is not important in the low-gradient limit. 

In the present section we discuss modifications of the 
picture drawn above that arise due to the electron- 
electron interaction. From the perspective of the 
"fermionic" solution of Sec. Ill Bl one obvious consequence 
of the interaction is the appearance of energy relaxation 
leading to local thermalization of the distribution func- 
tion. This thermalization will eventually wash out all the 
oscillating features of the density. However, the corre- 
sponding time will be very large, since the lifetime of elec- 
tronic excitations in an interacting ID system scales as a 
high power of the mass m (inverse curvature of the spec- 
trum near the Fermi points), or equivalently, of the Fermi 
momentum pf — tuvf, see Ref. 8 for a review. Specifi- 
cally, at zero temperature the lifetime of a quasiparticle 
with momentum k due to a long-range (smooth on the 
scale Xp) electron-electron interaction V(r) is given byS^ 

k F ) 4 



1 



[V (Vo - V k - kp )] 



(k 



(58) 



where V q is the Fourier transform of V(r). We will be 
particularly interested below in the case of power-law de- 
caying interactions 

nr) = -isri, (59) 



mln 



for which Vq — V„ cx q° 



Here, 1 < a < 3, and the 



length Iq parameterizing the strength of the interaction 
is the Bohr radius for the potential V(r). Estimating 
now the relevant momentum k as k — kp ~ Ap we find 
the inelastic decay rate^ 

1 - 1 r aoPoo) 2Q - 2 (w) 2 + 2 ". 



m 7 v F lQ 



(60) 



If the interaction falls off faster than 1/r 3 , one has 
Vo — V q oc q 2 ; the corresponding result can be obtained by 
setting a — 3 in Eq. (|6"0")) . On the other hand, the charac- 
teristic time scale for the density ripples is the shock time 
t c ~ mAs/Ap. Assuming moderate interaction strength 
lo 



1/ Poo we find 



N 



Ap 

Poo 



(61) 



We see that in the limit of small Ap/ ' p^ <C 1 the char- 
acteristic time t p of inelastic decay given by Eq. (IBTJl) is 
much larger than the shock time t c . In other words, the 
relaxation effects remain negligibly small at times much 
larger than t c . In view of this, in the rest of the pa- 
per we neglect the influence of inelastic relaxation on the 
dynamics and focus on other interaction-induced effects 
that strongly affect the development of density oscilla- 
tions. 
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A. Finite-range interaction 

Let us first briefly discuss the influence of finite-range 
interaction on the dynamics of fermions. We parametrize 
the interaction potential at low momenta by the scatter- 
ing length Iq and the effective interaction radius l{ n t 



V q = -±-(l-q% 



zr-^» (62) 



1 



where l\ = l 2 nt /l$. Correspondingly, the interaction- 
induced correction to the Hamiltonian takes the form 

(63) 



H, 



int 



P T -h (d xP f 



Within the hydrodynamic description the fermionic mass 
m manifests itself only via the time scale t c and we have 
set m to unity (cf. the case of free fermions, Sec. |TTJ) . We 
see that the zero momentum component of interaction 
gives rise to an additional p 2 term in the Hamiltonian. 
The only effect of this correction is the renormalization of 
Fermi velocity. The (j 2 -part of the potential V(q) renor- 
malizes the p 2 term in the free Hamiltonian. The re- 
sulting term may compete with the last term (the one 
containing the Hilbert transform) of Eq. (|57l) in govern- 
ing the dispersive regularization of the shock dynamics. 
If the interaction range is not too long, l\ <C 1/Ap (this is 
in particular the case for a short-range interaction with 
h -Vf), the interaction- induced p 2 term can be dis- 
carded, and the dynamics will be the same as in the free- 
fermion case. In the opposite limit of a very-long-range 
interaction, 1% ^> 1/Ap, it is the the interaction-induced 
term that will control the dispersive regularization. Con- 
sequently, the period of oscillations will not be given any 
more by the free-fermion result (|56j) but rather will have 
a form of Eq. (|3"3")l with A^ replaced by 1%, which yields 



Sx ~ v / ^i/A / 9> 1/Ap. 



(64) 



Corresponding equations and their solutions are dis- 
cussed in Appendix [3] the only difference is that the 
dispersive term is now enhanced by a factor ~ Pooh- 
Similarly to what we will see below for power-law inter- 
actions ( Sec. IIV SI and ITVCl) . the character of resulting 
oscillations will now depend on the sign of the initial 
pulse. Let us assume that the interaction is repulsive. 
Then for an initial density dip the oscillation will have a 
shape similar to those of free fermions [but with a larger 
period according to Eq. (|64[)] . On the other hand, for 
an initial hump, the perturbation will decompose in a se- 
quence of well-separated solitons, cf. Fig. [3] below. The 
particle number q carried by each soliton is obtained from 
Eq. (|A18[) by a replacement l/poo — > h, which results in 
q ~ v/ZiAp > 1. 



B. Calogero model 

The Calogero-Sutherland (CS) mode l 30 ' 31 is a remark- 
able example of the quantum integrable model. It 



appears in various branches of physics, such as spin 
chains, disordered metals and fractional quantum Hall 
edge s 12 ' 32 " —. 

In the CS model the particles interact via an inverse- 
square potential 



V(x) 



A(A-l) 



mi' 



(65) 



Here A is the dimensionless interaction strength and m 
is particle mass. We will confine ourselves to the case of 
strong repulsion, A ^> 1. 

Being interested in the hydrodynamic description of 
the CS mode l 11 ' 36 " —, we have to rewrite the CS Hamil- 
tonian in terms of the particle density. While the free 
part of the Hamiltonian after bosonization turns into the 
cubic Hamiltonian (|21|) . we need also a regularized ex- 
pression for the interaction term 



H-mt = - I dxdx'V(x 



x')p{x)p{x') . 



(66) 



The necessity of the regularization arises due to singu- 
larity of V(x) at x = and the corresponding ultraviolet 
divergence of the interaction at zero momentum. Taking 
the inverse particle density as the natural ultraviolet cut- 
off in the problem, we can rewrite the interaction term 
as 



int. 



A(A - 1) / dxp' 



A(A-l) 



dxpA-ip . (67) 



The operator A2 was defined in Eq. (|35|) . The precise 
coefficient in front of the cubic term entering H lrA is out 
of control within this estimate. Also, terms with higher 
gradients of the density may appear upon accurate reg- 
ularization of the model. A more rigorous treatment of 
the CS model^^i - leads to a slight modification of (JBTJ) 
and results in 



H 



dx 



pv 



tt 2 A 2 p 3 A(A - 1) 



(A 



pA 2 p 



1) 2 P 21 



(68) 



The characteristic feature of the Calogero model is the 
scaling of the interaction with distance which coincides 
exactly with that of the kinetic energy. At large coupling 
constant, A > 1, the potential energy dominates over 
the kinetic energy at all scales and drives the model to- 
wards the semiclassical limit. Indeed, rescaling by A the 
density and the space-time coordinates and switching to 
Lagrangian formalism, one finds the action correspond- 
ing to Eq. (f68|) : 



A / dxdt 



1 



b 2 

XX 



7T~ 
"if 



+ 



1 



(69) 

Here d x (f> = p. The large factor A in justifies now 
the semiclassical approach. 
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FIG. 9: Density profile for the Calogero model with strong 
repulsive interaction after the shock, t — 5t c . Initial pertur- 
bation was a density hump (top panel) and density dip (bot- 
tom panel). In the case of a density hump, the initial gaus- 
sian density perturbation decays into solitons carrying exactly 
one electron each. On the other hand, for a downward density 
pulse in the initial state, nearly sinusoidal oscillations develop 
in the shock region, so that the density evolution is similar to 
that of free fermions. Note that for a density dip the shock 
occurs on the rear side of the pulse. 



analytically via the solution of modulation equations dis- 
cussed previously in Sec. lIII Cl One finds (see AppendixlO 
for details) that close to the front edge of the train the 
height and width of the solitons (which are of Lorentzian 
shape) are given by 



Sp = 2Ap , Sx — 



1 



2irAp 



(70) 



The density evolution is very much different (and much 
more similar to that of free fermions) for the case of an 
initial density dip. The corresponding data are shown in 
the bottom panel of Fig. [9j they are fully analogous to 
the previous results of Fig. [5] (Note that for a density 
dip the shock occurs on the rear side of the pulse.) We 
see that a nearly sinusoidal dispersive wave is formed in 
the shock region. 

Such a dramatic difference in the behavior of particle- 
like and hole-like pulses has a simple qualitative expla- 
nation. A soliton can be formed when the effects of non- 
linearity and dispersion on the velocity counteract; their 
balance yields a soliton that moves preserving its shape. 
In the case of the Lagrangian (|69|) [which is equivalent, 
up to an overall factor A, to L/, of Eq. ([37]) ] the disper- 
sive term reduces the velocity, see Eq. ([31)]) with the lower 
sign corresponding to L^. Therefore, solitons can form if 
the non-linear term will enhance the velocity. This is the 
case when Ap is positive, i.e. for a density hump. 



C. Coulomb and other slowly decaying interactions 

Let us now turn to interactions decaying slower than 
the inverse distance squared (for definiteness, we will as- 
sume a repulsive interaction), 



Va(r)= 2 



1 1 



(71) 



Note that the Lagrangian in Eq. ([69[) is precisely the 
"hole" (not particle!) Lagrangian Lh we encountered in 
our discussion of loop corrections to hydrodynamics of 
free fermions [see Eq. ([37]) of Sec. lIII C] . The correspond- 
ing equations of motion are the Euler equations pip with 
the enthalpy given by (|39|) except for additional minus 
sign in front of the non-local term. 

This change of sign has a dramatic effect on the den- 
sity evolution in the system after the shock, as illustrated 
in Fig. [9] (top panel) where we plot the fermionic density 
at t ~ 5t c for the same initial density hump as was used 
previously. In the shock region, instead of the disper- 
sive wave seen in Fig. [8l one observes the formation of a 
solitonic train. Thus, at late stages of the evolution the 
initial hump decays into well separated solitons. Each 
of the solitons carries exactly one particle. The quanti- 
zation of solitonic charge, which is equal to unity, is a 
distinct feature of the strongly repulsive Calogero model. 

The solitonic train in the shock region can be studied 



with an exponent a satisfying 1 < a < 2. The case a = 1 
corresponds to the Coulomb interaction and is the most 
relevant from the experimental point view. Throughout 
this section we will assume the interaction to be weak in 
the sense that the parameter r s = Xf/Iq is small, r s <l. 
(There is no problem in analyzing the strong interaction 
regime, r s > 1, in a similar way, and we expect a quali- 
tatively similar behavior.) For a = 1 we will also assume 
that the Coulomb interaction is screened at a sufficiently 
large distance d^> Iq. 

The reasoning of the previous section which led us to 
Eq. ([57)1 is easy to generalize for the present case with 
the result 



/./;,„ „ " — / dx [ p a+1 - l - P A aP ) . (72) 



A Q (q) = -2r[l- a] sin^r 1 



a > 1 , (73) 
Ai(g)=loggd, gd»l. (74) 
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FIG. 10: Density profile after the shock for the case of Coulomb interaction. The initial density perturbation is positive and 
coincides with the one used in Sec. [II] (up to the constant equilibrium density poo). The Bohr radius lo, the screening length 
d, and the snapshot time t are indicated in each graph. For convenience of presentation, a constant coordinate shift has been 
made: x — > x — 20000 for the left (5t c ) plots and x — > x — 30000 for the right (7t c ) plots. For the upper two and the middle 
two graphs, the half-maximum width of the main peak is approximately given by d, as shown in the right plots. For the two 
bottom graphs the width is smaller than d by a factor ~ 1.5. This deviation from scaling with decreasing d is probably related 
to the fact that the condition II 7811 becomes less well satisfied. 



In the Coulomb case the p a+1 term in -ffmt.a should be 
replaced by (p 2 log dp)/ lo. For a > 1, the precise nu- 
merical coefficient in front of the p a+1 term can not be 
found within this reasoning. On the other hand, this 
term is small in parameter r s compared to the cubic term 
in the Hamiltonian of free fermions which provides the 
dominant nonlinearity. The only effect of the non-linear 
contribution iJj nt is a renormalization of Fermi velocity 



(small at r s -C 1) and we can omit it. Combining the 
free bosonised Hamiltonian with the relevant (dispersive) 
part of the interaction correction, we find 



H 



dx 



pv 



— - 2 PAaP 



(75) 



As we have found previously, the loop contribution to 
the equations of motion can be modeled by a correction 
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to the Hamiltonian of the form 
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FIG. 11: Density profile after the shock for the case of 
Coulomb interaction. The initial density perturbation differs 
from the gaussian pulse of Sec.[ll]by the change of sign. The 
legend indicates the Bohr radius io and the screening length d. 
For convenience of presentation, a constant coordinate shift 
x — > x — 30000 has been made. In contrast to the case of a 
positive (upward) density pulse (Fig. I10[) . a nearly sinusoidal 
oscillatory behavior is observed in the shock region. 
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FIG. 12: Density profile after the shock (time t ~ 7t c ) for 
electrons interacting via r~ 3 ^ 2 potential. The initial density 
perturbation is positive and coincides with the one used in 
Sec. [IT] (up to the equilibrium density poo)- The Bohr radius 
lo is indicated in the plots. For convenience of presentation, 
a constant coordinate shift x —¥ x — 120000 has been made. 



SH ioop ~ / dxpA 2 p 



(76) 



Comparing the interaction-induced dispersive contribu- 
tion [last term in Eq. (|75l) ] to Eq. (|76|) . we see that the 
interaction controls the dispersive effects at scales larger 
than Iq. 

The characteristic scale developed by the density per- 
turbation after the shock results from the trade-off be- 
tween non-linearity and dispersion. For a > 1 the corre- 
sponding estimate yields the scale 



5x ~ Iq 



1 



(^oAp) 



P 



(3 = 



1 



(77) 



We see that Sx is indeed much larger than Iq provided 
that 

l Ap < 1 . (78) 
In the Coulomb case under the same assumptions we get 



5x oj d 3> lo ■ 



(79) 



Thus, under the assumption (1781) the physics of density 
oscillations is indeed dominated by scales much larger 
than Iq, so that the neglect of loop corrections is justified. 

Experience gained in the analysis of the Calogero 
model and free fermions allows us to predict qualitative 
features of the density evolution, most prominently, its 
dependence on the sign of the density perturbation and 
the sign of the interaction. Specifically, we expect that for 
repulsive interaction and positive perturbation a train of 
solitary waves should emerge in front of the pulse. Each 
solitary wave is expected to carry 5xAp 3> 1 particles. 
On the other hand, a downward density perturbation (for 
the same repulsive interaction) will lead to formation of 
a nearly sinusoidal dispersive wave in the shock region. 
The change of the sign of interaction will result in the 
interchange of these two types of behavior. 

To support the qualitative analysis presented above, 
we have performed numerical simulations of the hydro- 
dynamic equations dictated by the Hamiltonian (1751) . Let 
us discuss the Coulomb case first. Figure QJJ shows the 
density perturbation for times t = 5i c and t = 7t c for 
electrons interacting via a Coulomb potential. The initial 
density hump was Gaussian with the same parameters as 
in the previous sections except for the equilibrium den- 
sity Poo which was taken larger to ensure that Iq > Xp^. 
Values of the parameters lo and d as well as of the time 
t are indicated in each of the plots. We clearly observe 
formation of a solitary wave and beginning of the for- 
mation of a second one (better pronounced for smaller 
d). According to the estimates presented above, in the 
Coulomb case the characteristic scale of the density oscil- 
lations emerging after the shock should be given simply 
by the screening length d. This is indeed confirmed by 
our numerics. In particular, the half-maximum width of 
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the main peak in the plots with d = 100 and d = 75 is 
equal to d. The two bottom plots (with the smallest d 
equal to 50) demonstrate some deviation from this scal- 
ing. This is possibly related to the fact that the condition 
(|78p becomes less well satisfied in view of increasing am- 
plitude Ap of the peak. 

Figure [TT] illustrates the change of the density behavior 
upon the change in of the sign of the density perturba- 
tion. As expected, for negative Ap we observe onset of 
nearly sinusoidal oscillations with a period ~ d in the 
shock region. 

We have also performed numerical study of fermions 
interacting via the intermediate potential V%i2(x). The 
results are exemplified in Fig. [T2l We observe that the 
density develops a solitary wave, similarly to the case of 
Coulomb interaction Fig. [TUJ The scaling of the width 
of the soliton agrees well with our above estimate of the 
characteristic scale for a = 3/2, 



if we use for Ap the actual amplitude of the peak. (While 
in the Calogero 1/r 2 case the soliton amplitude Ap is 
determined by that of the initial pulse, this is no more 
true for a < 2.) Note that the parameter loAp remains 
sufficiently small (~ 0.04 for the upper plot and 0.2 for 
the lower plot), so that the neglect of loop corrections is 
reasonably well justified. 

The analytical arguments and numerical data pre- 
sented above unambiguously show that a sufficiently 
strong and sufficiently long-ranged interaction dominates 
over quantum corrections in controlling the dispersive 
effects. In this limit the Hamiltonian (|75|) and corre- 
sponding hydrodynamic equations provide a controlled 
description of the non-equilibrium dynamics in a quan- 
tum many-body system. 

V. SUMMARY AND OUTLOOK 

In this article, we have explored the evolution of a den- 
sity pulse in a ID fermionic fluid. Our focus was on the 
regime of a wave "overturn" (population inversion) that 
is induced by spectral curvature. We showed that be- 
yond the corresponding time the density profile develops 
strong oscillations with a period much larger than the 
Fermi wave length and performed a detailed analysis of 
these oscillations. We have considered the case of free 
fermions as well as various interacting models, including 
a finite-range interaction, CS model, generic power-law 
interaction, and screened Coulomb interaction. Our key 
results can be summarized as follows: 

1. For the case of free fermions we have studied the 
problem by means of direct quantum simulations. 
Further, we have obtained analytical solution us- 
ing the phase space representation and the Wigner 
function. The Wigner function of the initial state 



exhibits oscillations in the phase space (as a func- 
tion of the momentum). When the initial perturba- 
tion is allowed to propagate (i.e. after the quench) , 
the curvature of single particle spectrum leads to 
the formation of inverse population of electrons at 
times t > t c . In this regime, the oscillations of 
Wigner function in phase space induce real space 
density oscillations, with each "ripple" containing 
a fraction of an electron. The characteristic period 
of these oscillations is controlled by the amplitude 
Ap of perturbation and is independent of the equi- 
librium density poo (or, equivalently, of the wave 
length Xp). 

2. We have also addressed the free-fermion problem 
using a hydrodynamic approach. The semiclassic 
equation of motion leads to formation of a shock in 
the regime where the inverse population of fermions 
in the momentum space is generated. This shock is 
regularized by gradient corrections to the Hamilto- 
nian and by quantum fluctuations. We show that 
for free fermions the latter effect is more impor- 
tant. We model the quantum correction by includ- 
ing in the theory a dispersive term corresponding 
to a particle or hole branch of the fermionic spec- 
trum, depending on the sign of the initial perturba- 
tion. This yields two different hydrodynamic theo- 
ries (with a difference in the sign of the dispersive 
term) for upward and downward density pulses. 

We show that this approach correctly captures 
the period of shock-induced density oscillation but 
overestimates their amplitude. In the hydrody- 
namic language, the formation of oscillations is 
caused by an interplay of the non-linearity (caused 
by the spectrum curvature) and the dispersion 
(dominated, in the case of free fermions, by quan- 
tum corrections carrying information about the 
spectral curvature). 

3. The electron interaction leads to additional disper- 
sive terms in the hydrodynamic equations. For in- 
teraction that decays with the distance r slower 
than 1/r 2 , such terms dominate the long-distance 
(small momentum) behavior, and quantum correc- 
tion can be neglected. In this case the applicability 
of semiclassical hydrodynamic equations becomes 
fully justified. The case of CS model (1/r 2 inter- 
action) is marginal; the interaction-induced disper- 
sive term is dominant (and thus the semiclassical 
hydrodynamic approach is fully controlled) if the 
interaction is strong, A> 1. 

For the case of a finite-range interaction the domi- 
nant dispersive term is provided by the interaction 
only if the interaction radius is very big; otherwise, 
the free-fermion results apply. 

4. In the situations when the interaction controls the 
dispersive effects (and thus the semiclassical hy- 
drodynamic approach is fully under control), the 
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impact of interaction depends on its sign and 
the sign of the density perturbation. Specifically, 
for a repulsive interaction and a density dip (as 
well as for an attractive interaction and a density 
hump), we observe formation of nearly sinusoidal 
oscillatory structure similar to the free fermions 
case. Quantitative characteristics of the oscillations 
(wave length and a number of particle in each "rip- 
ple" ) are however in general parametrically differ- 
ent compared to the free-fcrmion model. 

On the other hand, for a repulsive interaction and 
a density hump (as well as for an attracting inter- 
action and a density dip), the interaction leads to 
the formation of a train of solitary waves. In gen- 
eral, the charge (particle number) carried by each 
soliton is non-universal (depends on the type and 
the strength of the interaction, and on the ampli- 
tude of the perturbation). A notable exception is 
the CS model with A ^> 1, when the solitons carry 
a unit charge. 

We hope that our predictions can be verified experi- 
mentally. There is a number of electronic realizations of 
ID fermionic systems, including carbon nanotubes, semi- 
conductor and metallic nanowires, as well as quantum 
Hall and topological insulator (quantum spin Hall) edges. 
For these electronic liquids, a model with Coulomb inter- 
action is expected to be applicable (except if special ef- 
forts are made to strongly screen it) . An alternative phys- 
ical realization is provided by systems of cold fermionic 
atoms. This is probably the most natural experimental 
realization of the models of free fermions and of finite- 
range interaction. 

Before closing, we list some of directions of further the- 
oretical research opened by the present paper; a work in 
some of these directions is currently underway. 

1. An interesting question is whether it is possible to 
formulate a more general classical hydrodynamic 
theory for free fermions that would controllably 
capture evolution of a generic density perturba- 
tion, including both upward and downward density 
pulses. Such a theory can be useful from the funda- 
mental point of view, as well as for the problem in 
which quantum corrections and interaction effects 
are comparable. 

2. For models with power-law interaction other than 
CS model (including the experimentally most rel- 
evant case of the Coulomb interaction), it is im- 
portant to complete analytical investigation of the 
emerging oscillations and solitary waves and to ex- 
plore the integrability of these theories. 

3. An important task is to perform ab initio calcula- 
tions for many-body quantum interacting system. 
The results should allow one to verify our above 
predictions (obtained in the framework of the hy- 
drodynamic theory) and to explore the interplay of 



quantum corrections and interaction (e.g. in the 
model with a finite-range interaction). 

4. Our results on evolution of a density perturba- 
tion should be also relevant to strongly repulsive 
ID bosonic problems, in particular, in view of the 
equivalence between the Tonks-Girardeau gas and 
free fermions. It would be very interesting to study 
the crossover from the quasi-condensate regime 
characteristic for weakly interacting boson a 16 ' 42 to 
the Fermi-like behavior for strong repulsion. On 
the experimental side, such a setup can be realized 
in the framework of cold bosonic atoms. 
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Appendix A: Hydrodynamics denned by Eqs. (|31|). 
([32]): Solitons and periodic solutions 

In this Appendix we analyze properties of the classical 
hydrodynamics defined by Eqs. (pTTj) . (|3"2"|) . Such a theory 
arises if we take into account dispersive terms generated 
by phenomenological Haldane's formalism, see Sec. MI Al 
but neglect the quantum loop correction. 

As explained in the main text, this turns out to be 
not a correct description of free fermions (since the loop 
corrections generate parametrically more important dis- 
persive terms). Nevertheless, the analysis of this theory 
is quite illuminating, and we present it in this Appendix. 
Furthermore such a theory arises in a fully controllable 
way in a model of finite-range interaction with a suffi- 
ciently large interaction radius, see Sec. IIV A[ 

Let us focus on a traveling wave excitations 

p(x,t) =p(x-Vt), (Al) 
v(x, t) = v(x - Vt) . (A2) 

Substituting Eq. (fATj) into Eq. (|3T]), one finds 

-Vd x p + d x (pv) = 0, (A3) 

d x y - Vv + \ 2 + w^j = 0. (A4) 

We now analyze some simple excitations described by 
Eqs. dMI, dSS- 
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Single soliton 



The Euler equation, Eq. (IA4[) . can thus be rewritten as 



We start with a solitonic wave. In this case the exci- 
tation of the density and velocity fields are confined to a 
finite region in space, and the continuity equation (|A3[) 
yields 



v = V 



P~ Pa 
P 



(A5) 



Substituting Eq. 
tains 



V 2 P 2 



P 



5J) into Euler equation (|A4I) . one ob- 
1 



tt 2 P 2 



- 2 pI 



:d x 



(A6) 



Since both sides of the equation are full derivative with 
respect to x, the order of this equation can be easily 
reduced, yielding 



P.,; 



^ 2 ploP 2 + ^Ep 



(A7) 



Here E = p oo T/ 2 /2 + 7r 2 y0^ o /6 is a constant of integration. 
Defining £ = p — prx, , one obtains 



Wa + bii + ce 



2x. 



(A8) 



where a = it 2 p 2 ^ — V 2 , b — 4tt 2 Poo/3, c = 7r 2 /3. Perform- 
ing the integral over £, we find the solitonic solution 



(z - b) 2 - Aac 



z = 2a + ^" c -2y^|x| 



(A9) 
(A10) 



where £o < is the largest (smallest by absolute value) 
root of a + b£ + c£ 2 . We note that soliton propagate 
with the velocity smaller than the velocity of sound (V < 
tt/Ooo), i.e. is a hole-like excitation from the fermionic 
point of view. The charge of a soliton 



y/3. 2 

Z7T 



3(1 — V 2 ) 



V 



V 



2+^3(l-V 2 ) npoc 
is less than unity and is not quantized. 

2. Periodic wave 



(All) 



Periodic solutions of Eqs. (| A3 j) ■ (|A4j) can be conve- 
niently parametrized by A4 > A3 > A2 > Xv 

4?r 2 

{d xP f = —(P - Xi)(p - A 2 )(p - Aa)(p - A 4 ) , (A12) 

where parameters A satisfy the constraint Ai + A2 + A3 + 
A4 = 0, and are related to the velocity of the wave ac- 
cording to 

V 2 = y (A! + A 2 + A 2 + A 2 A 3 + A 2 A 4 + A3A4) . (A13) 



dp 



Air 2 



dx . 



V / (A4-P)(A 3 - P )(p-A 2 )(p-A 1 ) 
Integrating this equation, we find 

Ai(A 4 - A 2 )[l - dn 2 (y; k)] - A 2 (A 4 - A x ) 



(A14) 



p(x) = 



(A 4 -A 2 )[l-dn 2 ( 2 /;fc)]-(A 4 -A 1 ) 



(A15) 



Here dn(y; k) is Jacobi elliptic function, 
y — x\J (A 4 — A 2 )(A 3 — Ai)7r 2 /3, and the elliptic modulus 



k is given by 



(A 3 -A 2 )(A 4 -A!) 
(A 3 - A1XA4-A2) 



(A16) 



The density oscillates within the interval A 2 < p < A 3 . 
The limit A 4 < A 2 < A 3 — > A 4 corresponds to trains 
of well separated dips (which are nothing but solitons 
considered above), where the size of each dip d is much 
shorter then the distance between the neighboring dips 
L. In this regime the width of the dip can be related to 
the density amplitude Sp as 



1 



vSppo 



(A17) 



The number of particles carried by a dip can be estimated 
as 



q = d5p - y/Sp/poo < 1 



(A18) 



The limit Ai < A 2 — > A 3 < A 4 describes a 
small-amplitude periodic wave. The amplitude of the 
w ave is s A 3 — A 2 , and its wave length is L = 
y/3/47t 2 (A 4 — A 3 )(A 2 — Ai). Thus the number of elec- 
trons carried by each "ripple" (period) of the wave is 



A3 — A 2 



v/(A 4 -A 3 )(A 2 




(A19) 



It is worth mentioning that the theory considered in 
the present Appendix bears a close connection with the 
KdV equation of the classical hydrodynamics. This is be- 
cause the regularizing term in the hydrodynamic equa- 
tions has a similar (third-derivative) structure in both 
cases. 



Appendix B: Modulation equations for 
hydrodynamics denned by Eqs. (f3T|). (f39]) 

In this section we address the issue of modulation equa- 
tions for the hydrodynamic system (pTTj) , (|59")l . Within the 
framework of Whitman modulation theory, one promotes 
the single-phase (periodic) wave (|40)) to an Anzats 



<i> = e(x,t) + ®(e(x,t),x,t) 



(Bl) 
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FIG. 13: Comparison of predictions of the modulation theory 
(dashed line) and numerical simulations (full line) for the hy- 
drodynamic theory (|3ip . (|39p . A perfect agreement between 
the analytic and numerical results is observed. 



and identifies the parameters of the single-phase wave 
with the derivatives of the phases 8 and 8: 



Pa = d x 8 , 
k = d x 8, 



7 : 

UJ 



-dtO 
-Stt 



(B2) 
(B3) 



Obviously, parameters defined in this way satisfy the con- 
tinuity equations 



dtPo + dxl = , d t k + d x uj = . 



(B4) 



To derive the modulation equations, we substitute the 
single- phase wave (gDJ), (|1T|). (|32"]l. (|15|1 into the La- 
grangian L p , neglect derivatives of the modulation pa- 
rameters, and average the result over a period of oscilla- 
tions. We get 



(L p )=--tt 2 pI + ^- 
6 1 2pi 



fc 3 
48^ 



UJ 

ink 



7^ 
27rpi 

2 



8tt 2 pi 
7fc few 



(B5) 



Here, a = sign(fc7 — ujpo) and pi is given by (|43j) . We 
now vary the averaged Lagrangian (|B5|) with respect to 
phases 8 and 8, keeping in mind the relations (|B2[) . (|B3I) . 
This yields 



d x 



d(L P ) 

dk 
d(L P ) 

op 



•a 



duj 
d(L p ) 



= 0. 
= 0. 



(B6) 
(B7) 



These two equations, together with the continuity equa- 
tions, constitute four equations for the four unknown 
parameters. Writing them explicitly and performing a 
change of variables according to Eqs. (|4*5|) . (fl9"l) . 

([SO]) , we arrive at Eq. Pol - 



Figure [T3l demonstrates the density in the shock region 
predicted by modulation theory together with the result 
of numerical simulations. We observe a perfect agree- 
ment between the analytical and numerical results. 



Appendix C: Modulation theory and soliton trains 
in Calogero model 

In Appendix [Bl we have discussed in detail the modu- 
lation theory of hydrodynamic equations (|3Tt , (|39t gen- 
erated by the "particle" Lagrangian L p , Eq. (p?Tl) . Let 
us now briefly address the modulation equations for the 
theory defined by the "hole" Lagrangian Lh and their 
solution for the upward density perturbation in the ini- 
tial state. This issue is relevant for the description of the 
solitonic train emerging from the positive density pertur- 
bation in the repulsive Calogero fluid (see Sec. HVBj) . 

The starting point for the modulation theory is a 
single-phase periodic wave. Its form can be obtained 
from Eqs fH]), (Q, (03) via the replacement k ->• -fc 
and uj — » — uj. The modulation equations can now be 
derived exactly in the same way as in Appendix [B] The 
result reads 



Po = 



k = u 3 -u 2 , 

-UQ + Ul - U 2 + u 3 

2^ ' 



-Un + Ui — U% + ItS 



47T 



with the Riemann invariants satisfying 
d t Ui + UidxUt = . 



(CI) 
(C2) 

(C3) 
(C4) 

(C5) 



Finally, applying the boundary conditions ([51"j) . (j52")l 
at the edges of shock region one finds the Riemann in- 

(i) 

variants in term of the branches of Fermi momentum p F 
(see Fig. EJ) 



P F 



i = 
u 



1,2,3, 

= -Poo ■ 



(C6) 
(C7) 



Here we have chosen labeling of Riemann invariants 
such that it exactly corresponds [see Eq. (|C6|) ] to our no- 
tations for Fermi momentum branches, as was also the 
case for the Lagrangian L p , Eq. (f53"| . Note that the re- 
lations between the parameters k, uj, p , 7 and the Rie- 
mann invariants differ in the two cases by cyclic permu- 
tation of Mi, u 2 and it 3 , cf. Eqs. (gTJ), (gSJ), (p[9"]l. (|5D)l and 
Eqs. (JUT]), (rC2|, JC3]), (IC41) . 

Contrary to the case of "particle" Lagrangian L p con- 
sidered in the previous section Eqs. (IC1|) . (IC2[) . (IC3|) 
and (|C4[) predict that the wave vector vanishes at the 
leading edge xi and the density perturbation decays into 
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Lorentzian-shaped solitons 

, . 1 A 

p[X,t) = poo + 



A = 



7T A 2 + {X - Vtf 
(!) i 

P F + P00 



V 



(3) , (2) 

P F +P F 



(C8) 
(C9) 

(CIO) 



In the limit Ap <C p<x> this simplifies to 



p(x,t) = poo + 



2Ap 



l + 4Tr 2 Ap 2 (x-V F t) 2 



(Cll) 
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